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ABSTRACT 


In  this  report,  the  requirements  are  considered  for  a  two-dimensional 
flow  visualization  diagnostic  capable  of  nonintrusively  measuring  tempera¬ 
ture  and  species  concentration  in  a  turbulent,  atmospheric  pressure  flow 
tube.  The  need  for  and  the  history  of  flow  visualization  techniques  are 
briefly  reviewed.  A  new  approach  is  then  suggested,  based  on  optical  tomo¬ 
graphy,  that  would  have  several  advantages  over  methods  currently  under 
development.  An  up-to-date  reference  list  for  optical  tomography  is  in¬ 
cluded.  The  feasibility  studies  completed  to  date  are  presented  in  the 
bulk  of  the  report. 
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I.  INTRODUCTION 


A.  BACKGROUND 

Current  research  in  reacting  and  nonreacting  flows  is  motivated  by  a 
desire  to  understand  basic  processes  and  specific  practical  devices.  For 
nonreacting  flows,  applications  to  sonic  systems  include  aerodynamics  of 
fixed  and  rotating  airfoils,  gas  turbines,  and  turbulent  mixing  problems. 
High-pressure  nozzles,  rocket  propulsion  (chemical  and  electric),  and  shock 
waves  are  supersonic  nonreactive  flows  of  interest.  Reactive  flows  can  be 
classified  as  steady  state  or  turbulent.  Steady-state  flows  often  studied 
are  stabilized  premixed  and  diffusion  flames,  and  laminar  flow  reactors  and 
incinerators.  Examples  of  turbulent  reacting  flows,  however,  are  consider¬ 
ably  more  numerous  and  of  greater  practical  importance.  These  include 
nearly  all  coiranercial  combustors  used  for  power  generation  and  operated 
with  fossil  fuels;  coal  gasifiers;  internal  combustion  engines;  open  cycle 
magnetic  hydrodynamic  (MHD)  power  generators;  and  the  study  of  propulsion 
systems  (gas  turbine,  ramjet,  afterburner,  and  rocket).  To  understand 
turbulent  reacting  flows,  hydrodynamics  and  spatially  specific  chemical 
kinetics  must  be  understood  in  order  to  quantitatively  model  the  system. 

A  system  problem  of  immediate  interest  to  the  Air  Force  Space  Division 
(AFSD)  is  the  disposal  of  hydrazine  by  incineration.  Of  particular  inter¬ 
est  is  the  rate  of  competitive  chemical  reactions  involved  in  the  combus¬ 
tion  of  ^Hjp  which  is  very  sensitive  to  temperature.  This  temperature 
sensitivity  and  the  dependency  of  the  final  products  (specifically,  N0X 
pollutants)  on  the  flow  contours  in  the  mixing  process  require  a  detailed 
understanding  of  the  relationship  among  the  macroscopic  flow  parameters, 
temperature,  and  chemical  kinetics.  Computer  models  are  often  used  to 
further  this  understanding.  However,  these  models  require  the  solution  of 
the  conservation  equations  of  continuity  and  momentum,  and  a  scalar  (tem¬ 
perature  or  concentration)  for  a  high- turbulence  Reynolds  number  in  a  den¬ 
sity-weighted  (Favre-averaged)  form.  Solution  of  these  equations  often 
involves  making  empirical,  unjustifiable,  or  even  erroneous  assumptions. 
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For  example,  turbulence  chemistry  closure  requires  a  probability  density 
function  that  is  chosen  somewhat  arbitrarily;  the  turbulent  transport  pro¬ 
cess  is  often  cast  in  the  form  of  a  turbulent  eddy  diffusivity,  a  concept 
increasingly  under  attack;  and  finally,  the  modeling  of  dissipation  terms 
is  empirical,  at  best  (Ref.  1).  These  assumptions  require  that  a  consider¬ 
able  database  be  collected  on  systems  of  practical  importance  before  models 
are  implemented. 

The  history  of  the  diagnostic  techniques  used  to  probe  hostile  flows 
parallels  that  of  scientific  apparatus  in  general.  Originally,  intrusive 
probes  were  used  to  sample  gas  streams  (thermocouples,  Pitot  tubes,  and 
devices  to  extract  and  route  samples  to  on-line  chemical  analyzers).  These 
devices  were  slow,  limited  to  relatively  benign  environments,  and  more 
importantly,  often  perturbed  the  parameter  being  measured.  Optical  tech¬ 
niques  did  not  have  these  problems,  and  they  proliferated  with  the  advent 
of  the  laser.  Line-of-sight  techniques  (sodium  line  reversal,  absorption 
spectroscopy,  and  one-  and  two-color  laser  transmissometry  for  particle- 
ladened  flows)  were  the  first  optical  techniques  to  be  implemented.  As 
technology  improved,  spatially  precise  (point)  measurements  could  be  made 
using  laser-induced  fluorescence  (LIF),  spontaneous  Raman  scattering  (SRS), 
Rayleigh  scattering  (RS),  coherent  anti-Stokes  Raman  scattering  (CARS), 
laser  Doppler  velocimetry  (LDV) ,  and  Mie  scattering  (MS).  Further  sophis¬ 
tication  led  to  simultaneous  measurements  of  points  along  a  line  for  LIF, 
SRS,  and  RS  (and  marginally  for  CARS). 

If  photographic  techniques  are  neglected,  the  field  of  flow  visualiza¬ 
tion  began  when  technology  developed  enough  to  allow  the  two-dimensional 
implementations  of  some  of  the  preceding  processes.  The  relative  strengths 
of  these  processes  are  compared  in  Table  1 .  This  table  reveals  that  the 
principal  disadvantage  of  SRS  and  RS  is  their  lack  of  strength.  Several 
planar  Raman  scattering  experiments  have  been  attempted  (Refs.  2-5). 

Because  of  the  exceedingly  small  scattering  cross  sections  of  Raman  scat¬ 
tering,  these  experiments  usually  employ  multipass  cells  relaying  a  high- 
intensity  laser  beam  through  a  small  measurement  area  in  a  nonluminous 
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Table  1.  Cross  Sections  for  Various  Physical  Processes 


Physical  Process 

Differential 
Cross  Section 

p 

(cnr/sr) 

Particle  Scattering  (Mie) 

10  ym  diam 

IQ'7 

0.1  pm  diam 

10'13 

Atomic  Fluorescence 

Strong,  Visible 

CO 

1 

0 

1 

00 

1 

0 

Molecular  Fluorescence 

Simple  Molecules 

-=3- 

C\J 

1 

O 

1 

1 

0 

Rayleigh  Scattering 

N21  532  nm 

CO 

CM 

1 

O 

X 

t*— 

Rotational  Raman  Scattering 

N2,  all  Lines,  532  nm 

7  X  10“3° 

N2,  Strong  Line  (Including 
Population  Factor),  532  nm 

3.5  x  IO'31 

Vibrational  Raman  Scattering 

N2,  Stokes  Q-Branch,  532  nm 

3  x  10~31 
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environment.  Because  the  Raman  signal  can  be  overwhelmed  easily  by  emis¬ 
sion,  RS  is  limited  to  the  most  benign  laboratory  environments.  Planar  Mie 
scattering  (Refs.  6-9)  and  Rayleigh  scattering  (Ref.  10)  have  been  used  to 
measure  temperature  or  composite  density  fields  in  atmospheric  pressure 
flows.  MS  requires  that  the  flow  be  seeded  with  particles.  Inferring 
temperature  or  density  from  these  data  requires  many  assumptions  (e.g., 
constant  pressure,  homogeneous  optical  characteristics  of  the  gas  mixture, 
uniform  particle  seeding,  and  uniform  particle  size)  that  restrict  these 
techniques  to  carefully  controlled  environments.  Additionally,  neither  MS 
nor  RS  can  differentiate  among  chemical  species.  Despite  these  shortcom¬ 
ings,  interest  in  these  techniques  continues,  because  the  higher  signal 
strengths  allow  the  plane  of  radiation  to  be  translated  quickly  through  the 
flow,  permitting  an  electronic  camera  to  record  a  sequence  of  planar  images 
that  may  be  later  reconstructed  into  a  three-dimensional  view  (Refs.  8-10). 

Most  recent  work  in  flow  visualization  has  employed  planar  LIF  because 
of  its  strong  scattering  cross  section  and  because  it  is  species  speci¬ 
fic.  Hanson's  group  at  Stanford  has  been  most  actively  developing  this 
technique.  We  reference  only  a  few  of  his  recent  publications  (Refs.  11- 
15).  Hanson  has  recently  extended  this  technique  to  include  two-photon 
absorption  (Refs.  13,  14)  and  a  three-dimensional  image  of  flow  fields 
(Ref.  15). 

Despite  these  successes,  LIF  has  two  drawbacks.  First,  it  requires 
optical  access  outside  the  plane  being  studied.  Second,  and  more  serious, 
is  collisional  quenching  of  the  laser  excited  state.  Quenching  is  a  large 
effect  that  greatly  hinders  the  quantitative  interpretation  of  the  data, 
thus  limiting  the  use  of  this  technique  to  extremely  well-characterized 
environments.  (Quenching  coefficients  are  sensitive  to  temperature  and  to 
the  relative  concentrations  of  surrounding  species.)  In  poorly  character¬ 
ized  environments,  the  data  are  only  useful  for  study  of  turbulent  struc¬ 
tures. 
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B.  MULTI  ANGULAR  ABSORPTION  AND  TOMOGRAPHIC  RECONSTRUCTION 


Simple  optical  absorption  was  one  of  the  earliest  optical  diagnostic 
techniques  applied  to  gas  flows  but  has  received  scant  attention  in  recent 
years  because  of  its  lack  of  spatial  specificity.  Nevertheless,  this  tech¬ 
nique  has  several  highly  desirable  attributes:  its  cross  section  is  large 
— the  same  as  that  for  LIF;  quenching  is  not  a  problem;  its  data  are 
straightforward  to  interpret;  it  is  species  specific;  and  it  can  be  used  to 
measure  concentrations  and/or  temperature.  The  principle  of  optical 
absorption  is  based  on  Beer's  law 

L 

-ln(I/I  )  =  J  n.(T)q.(v)dJt 
o 

where  IQ  (I)  is  the  intensity  of  light  before  (after)  it  traverses  the 
sample,  ni(T)  is  the  number  density  of  species  i  at  temperature  T,  q.(v)  is 
the  absorption  cross  section  at  frequency  v,  and  L  is  the  path  length  over 
which  absorption  occurs.  For  low  concentrations,  this  expression  can  be 
linearized 

L 

(I  -  I  )/I  =  J  ni(T)q(v)di 

o 

The  dynamic  range  of  optical  absorption  can  be  expanded  by  viewing  two 
spectral  lines  simultaneously.  The  strengths  of  the  absorption  cross  sec¬ 
tions  of  these  two  lines  are  chosen  to  differ  by  an  order  of  magnitude  (see 
Section  V.B). 

The  principal  disadvantage  of  absorption  measurements — the  lack  of 
site  specificity — can  be  eliminated  by  making  a  series  of  measurements  at  a 
number  of  different  angles.  This  collection  of  measurements  then  can  be 
used  to  mathematically  construct  a  cross-sectional  view  of  the  flow  (Fig. 
la).  This  procedure  is  the  principle  underlying  the  well  known  computer¬ 
ized  axial  tomography  (CAT)  scanners  used  for  medical  diagnostics.  The 
essence  of  the  procedure  can  be  understood  with  the  aid  of  Fig.  1.  The 
problem  consists  of  recovering  a  two-dimensional  function  F(x,y)  from  a 
series  of  projections  P(r,8)  that  are  related  to  the  absorption  as 
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TRAVERSE 

TYPICAL  SCANNER 


.  1  (a)  Geometrical  and  Mathematical  Relationships  Between  Projection 

P(r,8)  and  Cross-Sectional  Rays  F(x,y)  in  Tomography;  (b)  Geom¬ 
etry  of  a  Typical  Scanning  Parallel  Beam  Tomographic  Instrument 
(from  Ref.  20). 


P(r,e)  =  J  F(r,s)ds 

—oo 

00 

=  f  N.q.ds  =  -£n( I/I  )|  n 
“  iu  o  'r.fl 

—  CO 

The  rays  traversing  the  probed  region. need  not  be  parallel.  A  fan  beam 
geometry  (Fig.  2)  is  often  used  in  commercial  medical  CAT  scanners  because 
it  usually  results  in  shorter  data  acquisition  times. 

Inversion  algorithms  to  obtain  F(x,y)  from  the  experimentally  measured 
P(r,0)  are  generally  classified  as  either  Fourier  transform  or  linear 
superposition.  In  the  former,  the  one-dimensional  Fourier  transform  of  the 
projection  is  taken  to  obtain  the  spatial  spectrum  P(u,e): 

00  '  CO  CO 

P(w,e)  =  f  P(r,e)e-iurdr  =  J  J  F(r,s)e~iwrdrds 
—  00  —  00  —  00 

This  spatial  spectrum  can  be  written  as  a  two-dimensional  Fourier  transform 
using  the  "central  slice  theorem" 

00  CD 

P(u>,8)  =  J  J  F(r,s)exp[-i(4>r  +  ns)]drds|  _ 

—  CD  —00  ^ 

This  integral  is  evaluated  along  the  line  n  =  0  in  the  two-dimensional  fre¬ 
quency  plane.  The  spatial  frequency  is  given  by  p  and  us  =  2-np .  A  con¬ 
tinuous  sweep  of  180°  (360°  for  fan  beam  geometry)  is  required  to  define 
P(w,6).  The  desired  function  F(x,y)  is  obtained  as  the  inverse  Fourier 
transform  of  P(w,e),  which,  in  turn,  is  obtained  experimentally  from 
P(r,e).  This  procedure  can  be  implemented  quickly  using  fast  Fourier 
transform  techniques  (Refs.  16,  17),  but  a  "back  projection"  algorithm  has 
been  developed  (Ref.  18)  that  is  considerably  faster,  allows  the  use  of 
filter  functions  for  smoothing,  and  is  stable  when  P(r,e)  is  inferred  from 
noisy  data. 
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Fig.  2.  Geometry  of  the  Time-Multiplexed  Fan  Beam  Tomography  Instrument 
Introduced  in  the  Text 
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Linear  superposition  techniques  are  based  upon  algorithms  that 
repeatedly  adjust  the  estimated  F(x,y)  map  until  convergence  is  attained. 

An  example  of  this  procedure  is  the  algebraic  reconstruction  technique 
(ART),  which  divides  the  (x,y)  plane  into  a  grid  of  elements  that  further 
divide  each  beam  (r,e)  into  a  set  of  i  segments  of  known  weight  si 

P(r,8)  =  l  s.(x,y)  F.(x,y) 
i 

Iterative  methods  have  been  used  to  solve  systems  of  200  x  200  equations 
for  this  application  (Ref.  19). 

Optical  tomography  was  first  suggested  for  the  study  of  reacting  flows 
in  1980  by  Goulard  and  Emmerman  (Ref.  20).  However,  given  the  importance 
of  the  problems  that  optical  tomography  can  potentially  solve,  experimental 
implementation  of  the  technique  has  been  slow.  Goulard  and  Ray  (Ref.  21), 
and  Ray  and  Semerjian  (Ref.  22)  have  described  experiments  undertaken  at 
NBS.  In  these  experiments,  temperature  and  concentration  fields  were  meas¬ 
ured  in  a  symmetric,  steady-state,  premixed  methane/oxygen/diluent  flame 
using  OH  as  the  probed  species.  Narrow  frequency  light  produced  by  a  con¬ 
tinuous-wave  (cw)  ring  dye  laser  was  scanned  through  an  absorption  line  of 
OH  (see  Fig.  1b).  A  rotating  mirror  geometry  swept  the  beam  through  the 
burner,  creating  a  projection  in  1  ms  at  a  single  wavelength.  Scanning 
could  not  be  done  at  more  than  a  single  projection  angle;  therefore,  cylin¬ 
drical  symmetry  was  assumed  to  perform  the  reconstruction.  Absorption 
experiments  also  have  been  undertaken  at  Stanford  University  by  Byer  and 
his  students  (Refs.  23-25).  They  used  a  fan  beam  geometry  and  rotating 
mirror  to  create  reconstructions  of  cold  flowing  1 2  at  90  to  360  angles  on 
90  detectors  using  a  cw  Ar+  laser.  This  procedure  allowed  data  to  be 
acquired  for  reconstruction  of  concentrations  in  100  ms.  The  data  aquisi- 
tion  times  of  this  method  can  be  reduced  by  at  least  an  order  of  magni¬ 
tude.  Recently,  Synder  and  Hesselink  (Ref.  26)  demonstrated  a  novel  con¬ 
figuration  using  holographic  optical  elements  and  a  rotating  mirror.  In 
this  proof-of-principle  work,  only  one  angle  was  measured,  and  the  solid 
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sample  was  spun  to  create  the  multiple  projection  angles.  Variations  in 
the  index  of  refraction  can  be  used  to  measure  density,  and,  indirectly, 
temperature  fields  using  beam  deflection  tomography  (Refs.  27-30)  or  inter¬ 
ferometric  tomography  (Refs.  31-34).  All  the  work  cited  in  References  27- 
33  used  rotating  mechanical  elements  and  cw  lasers,  thus  restricting  the 
measurements  to  the  millisecond  time  scale.  This  time  scale  is  too  long  to 
freeze  turbulent  cells  in  unsteady  flows.  Recently,  however,  a  fast 
(300  us)  interferometric  tomographic  image  was  recorded  on  a  cylindrically 
shaped  piece  of  film  (Ref.  34).  This  demonstration  was  the  first  instan¬ 
taneous  optical  tomographic  reconstruction  of  a  nonstationary  fluid  flow. 
Several  implementations  of  optical  emission  tomography  have  recently  been 
demonstrated  (Refs.  35,  36).  Reviews  of  optical  tomography  for  applica¬ 
tions  to  fluid  flows  are  now  available  (Refs.  37,  38).  [Also,  see  the 
issue  of  Applied  Optics  devoted  to  computer  tomography  (Ref.  39).] 

C.  PROPOSED  INSTRUMENT  DESIGN 

A  novel  optical  instrument  for  tomographic  studies  of  turbulent,  re¬ 
acting  gas  streams  is  briefly  described  in  this  subsection.  This  instru¬ 
ment  differs  substantially  from  those  used  previously  in  experimental 
demonstrations  of  fluid  flow  tomography.  It  produces  higher  resolution 
(100  angles,  100  detectors /angle)  images;  has  faster  response,  thereby 
effectively  freezing  the  flow  (10-^  vs  10~^s);  uses  pulsed  rather  than  cw 
radiation;  and  has  no  moving  parts.  Furthermore,  the  instrument  is  easily 
scalable,  through  its  unique  time-multiplexing  design,  to  higher  resolu¬ 
tion,  greater  dynamic  range,  multispecies  detection,  or  even  to  a  config¬ 
uration  that  allows  the  temporal  development  of  turbulent  structure  to  be 
viewed  through  a  series  of  fast  snapshots. 

Figure  2  illustrates  the  geometry  of  the  proposed  instrument.  The 
center  circle  represents  the  cross  section  of  the  flow  to  be  investi¬ 
gated.  The  circumference  of  the  outer  circle  describes  a  locus  of  point 
sources,  each  of  which  consists  of  a  diverging  beam  of  light  described  by  a 
fan.  The  rays  from  each  source  (point  A  in  Fig.  2)  pass  through  the  probed 
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region  and  terminate  on  a  series  of  detectors  (d^-d^)  also  located  on  the 
outer  circle.  If  the  sources  are  well  separated  (points  B  and  C),  rays 
from  their  fans  will  not  overlap  on  the  opposite  side  of  the  circle,  and  a 
one-to-n  correspondence  will  exist  between  the  sources  and  detectors,  where 
n  is  the  number  of  detectors  in  each  fan.  Figure  2  depicts  a  geometry 
where  10  nonoverlapping  sources  could  be  placed  on  a  180°  semicircle.  If 
100  detectors  were  placed  in  each  fan,,  a  moderate  resolution  tomographic 
reconstruction  of  a  flow  could  be  performed. 

Now,  consider  placing  nine  additional  sources  between  sources  D  and  E, 
as  shown.  Further,  consider  that  these  sources  are  activated  sequentially 
so  that  source  t^  is  emitting  a  fan  of  radiation  at  time  t-|,  source  t2  is 
emitting  a  fan  of  radiation  at  time  t2>  etc.,  and  that  the  time  interval 
between  each  pulse  of  radiation  is  greater  than  the  time  any  source  is 
active.  Then,  a  detector  will  receive  the  radiation  from  only  one  source 
at  a  time.  A  train  of  these  radiation  pulses  can  be  created  from  a  single 
output  pulse  of  a  tunable  laser  by  relaying  the  pulse  back  and  forth  10 
times  between  two  retroref lectors,  laterally  displacing  it  at  a  fixed  dis¬ 
tance  on  each  pass.  A  partially  reflecting  mirror  inserted  at  45°  into 
this  ladder  of  parallel  paths  will  cause  a  portion  of  the  pulse  to  be 
directed  by  the  beam  splitter  through  a  lens  and  focused  into  a  bundle  of 
10  optical  fibers  on  the  first  pass  through  this  system.  These  10  fibers 
carry  their  simultaneous  pulses  to  the  10  locations  marked  t^  on  the  cir¬ 
cumference  of  the  outer  circle.  The  next  pass  through  the  beam  splitter 
creates  another  pulse  at  t2,  which  is  relayed  to  the  10  locations  on  the 
outer  circle  adjacent  to  t^.  In  this  way,  100  sources  are  created,  but 
only  10  widely  separated  sources  are  active  at  any  instant.  In  this  exam¬ 
ple,  each  source  illuminates  100  detectors,  creating  100  angle  x  100  detec¬ 
tor  resolution  with  only  1000  detectors.  Through  the  multiple  use  of  each 
detector,  a  larger  detector  size  is  possible  for  a  given  spatial  resolu¬ 
tion.  This  allows  the  acquisition  of  data  with  a  higher  signal-to-noise 
ratio. 
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Note  the  versatility  of  this  time-multiplexed  design.  By  the  intro¬ 
duction  of  a  second  sequence  of  pulses  into  the  fiber  at  a  different  wave¬ 
length,  a  second  spectral  line  can  be  probed,  allowing  both  the  temperature 
and  the  concentration  of  the  species  to  be  measured.  A  third  sequence  of 
pulses  allows  the  concentration  of  another  species  to  be  measured,  etc. 

Or,  the  second  sequence  of  pulses  can  be  used  to  obtain  a  second  "snapshot" 
of  the  same  species  at  a  later  time,  thus  tracing  the  time  evolution  of  the 
flow.  Multiple  snapshots  can  then  be  used  to  create  a  minimovie.  Current¬ 
ly,  an  integrated  circuit  is  available  that  allows  the  sampling  of  a  train 

4t 

of  640  pulses  at  time  intervals  of  10  ns  or  longer.  These  640  pulses  can 
enhance  resolution,  sample  multiple  species,  or  acquire  a  fast  series  of 
images . 

The  necessary  time  resolution  of  the  instrument  is  determined  by  the 
requirement  to  "freeze"  a  resolution  element.  The  resolution  element 
should  approach  the  smallest  space  scale  of  the  turbulence.  For  flow  vel¬ 
ocities  not  much  greater  than  100  m/sec,  the  smallest  space  scale,  or  eddy, 
will  be  approximately  1  mm  (Ref.  40).  If  an  element  is  considered  to  be 
frozen  when  it  does  not  move  more  than  10?  of  its  size,  all  data  must  be 
collected  in  10“^  sec.  As  the  velocity  decreases  and  the  scale  size 
increases,  the  time  resolution  of  the  instrument  can  be  relaxed  consider¬ 
ably. 

There  are  three  areas  of  preliminary  work  that  must  be  completed 
before  a  prototype  can  be  built:  (1)  analytic  studies  to  test  the  applica¬ 
tion  of  this  concept  to  a  transparent  tube  of  turbulent  gas;  (2)  software 
development  for  the  tomographic  reconstruction  algorithms;  and  (3)  labora¬ 
tory  studies  to  develop  components  of  a  prototype  system.  The  remainder  of 
this  report  presents  results  from  the  first  area  of  preliminary  work. 


* 

LeCroy,  monolithic  model  MVV200,  320/640-cell  analog  shift  register, 
LeCroy  Research  Systems  Corporation. 
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II.  GEOMETRICAL  CONSIDERATIONS 


A.  PLANAR  GEOMETRY 

As  noted  in  Section  I.B,  tomography  may  be  implemented  either  with  a 
parallel  beam  or  a  fan  beam  geometry.  The  concept  introduced  in  Section 
I.C  is  intrinsically  a  fan  beam  design.  Artifacts  introduced  from  the 
reconstruction  of  fan  beam  projections  limit  the  allowable  geometries,  as 
do  considerations  of  detector  size,  laser  speckle  noise,  minimum  size  of 
turbulence  cells,  etc.  Accordingly,  we  begin  analysis  by  deriving  a  number 
of  useful  relations  among  the  system  variables  for  a  planar  fan  beam  geom¬ 
etry  intercepting  a  cylindrical  flow. 

Consider  the  geometry  shown  in  Fig.  3.  A  fan  beam  of  full  angle  0 
emanates  from  point  B,  located  on  a  source-detector  circle  of  radius  R,  and 
passes  through  an  inner  concentric  circle  of  radius  r,  which  represents  the 
cross  section  of  the  flow  to  be  interrogated.  The  edges  of  the  fan  are 
tangent  to  the  flow  cross  section  and  strike  the  source-detector  ring  at 
points  C  and  D,  defining  an  arc  of  length  s.  The  arc  s  subtends  an  angle 
4>  with  the  common  center  point  A.  Noting 


2a  +  <t>  =  it 
28  +  0  =  it 


and  eliminating  a  and  8,  we  find 


*  =  20  (1) 
The  cosine  law  for  triangles  BCD  and  ACD  yields,  respectively 

c2  =  8a2  (1  -  cos0)  (2) 

c2  =  2R2  (1  -  cos*)  (3) 
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Fig.  3.  Geometry  Used  to  Calculate  the  Optimum  Planar  Fan  Beam  Configura¬ 
tion 
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Solving  Eqs.  (1-3)  for  4>(r,R)  and  noting  4>  =  s/R,  we  find  the  relationship 
between  the  arc  length  and  the  two  radii 

s  =  2R  cos” 1  ( 1  -  -r-x)  (4) 

R 

If  w  is  the  width  available  for  each  detector,  Nd  is  the  total  number  of 
detectors,  is  the  total  number  of  fibers,  Fm  is  the  multiplexing  factor, 
nd  is  the  number  of  detectors  per  fiber,  and  Na  is  the  number  of  arcs  per 
circumference,  then 


and 


2uR  =  N  ,w  =  N  s 
d  a 


N  .  =  n  ,Nr/F  =  N  n  . 
d  d  f  m  ad 


Combining  Eqs.  (4-6),  we  have 


R  =  {2/[l 
2uFr 


cos(tt/N  )]}^2  r  =  Fr 
a 


w  = 


N  n , 
a  d 


F  =  N./N 
m  i  a 


(5) 

(6) 


(7) 

(8) 
(9) 


where  F  is  the  ratio  of  the  diameter  of  the  source-detector  ring  to  the 
diameter  of  the  flow,  and  the  multiplexing  factor  Fm  is  the  number  of  times 
each  detector  is  used  in  one  reconstruction  data  set.  Multiplexing  allows 
for  a  lower  total  number  of  detectors  and  a  smaller  source-ring  diameter 
for  a  given  resolution.  For  large  Na,  the  factor  F  scales  as  Na, 
indicating  that  w  scales  as  r/nd,  independent  of  Na  and  the  multiplexing 
factor.  That  is,  for  a  given  resolution  (nd,Nf),  the  detector  size  is 
dictated  solely  by  the  flow  diameter.  This  constraint  can  be  relaxed  if 
the  source-ring  radius  is  not  equal  to  the  detector-ring  radius  (nonplanar 
geometry),  as  discussed  in  the  next  subsection. 
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Equations  (6) -(9)  can  be  used  to  select  an  appropriate  fan  beam  geom¬ 
etry  for  interrogating  a  cylindrical  flow,  subject  to  the  following  guide¬ 
lines.  Since  complexity  and  cost  of  the  proposed  tomographic  instrument 
are  directly  proportional  to  the  total  number  of  detectors,  but  the 
resolution  is  proportional  to  the  product  of  the  number  of  fibers  and  the 
number  of  detectors  per  fiber,  nd  should  be  as  large  as  possible,  and  Nd, 
as  small  as  possible.  This  implies  a  .large  factor  F  ,  which  in  turn 
requires  a  small  Nfl.  Na,  however,  cannot  be  made  too  small,  as  the 
practical  limit  for  the  fan  angle  lies  between  45  and  60°.  For  a  given 
flow  diameter,  a  compact  instrument  requires  a  small  source-detector  ring 
diameter  requiring  a  small  value  of  F.  There  are  two  other  reasons  for 
keeping  R  small.  First,  coherent  (laser)  radiation  passing  through  a 
multimode,  step- index  fiber  forms  a  speckle  pattern  on  the  detector  sur¬ 
faces.  This  random  pattern  is  a  source  of  noise;  therefore,  the  grain  size 
of  this  speckle  should  be  kept  small  compared  to  the  detector  dimensions. 
The  grain  size  increases  as  the  separation  between  the  fiber  and  detector 
increases.  Second,  beam  displacement  at  the  detector  surface,  resulting 
from  turbulence- induced  index  of  refraction  fluctuations  and  temperature 
gradients  within  the  flow,  increases  with  the  separation  between  the 
absorbing  element  and  detector.  Reconstruction  artifacts  also  place 
constraints  on  the  geometry.  For  example,  interpolation  error  is 
introduced  if  imaging  is  attempted  at  radii  less  than  0.5  R,  restricting  F 
to  values  greater  than  2  (Ref.  23).  Conversely,  the  detector  size  should 
be  kept  as  large  as  possible  (implying  a  large  R)  to  keep  the  preceding 
effects  small  compared  to  the  detector  dimensions.  Of  course,  the  larger 
the  R,  the  larger  nd  can  be  for  a  given  detector  width  w.  Thus,  the  ideal 
design  is  to  have  large  area  detectors  on  a  small  (but  not  too  small)  ring 
diameter . 

Finally,  for  the  commonly  employed  mathematical  expression  used  for 
the  fan  beam  geometry  in  the  convolution  back  projection  algorithm  [the 
Ramachandran-Lakshmunarayanan  (R-L)  kernel  (Ref.  18)],  spatial  resolution 
at  the  center  of  the  reconstructed  image  is  approximately  one-half  the 
detector  spacing  when  Nf  »  nd  (Ref.  23).  This  rule  of  thumb  allows  an 
approximate  number  of  pixels  to  be  assigned  to  a  given  fan  beam  geometry. 
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With  the  preceding  information,  the  fan  angle,  R/r,  the  diameter  of 
the  source-detector  ring,  the  detector  width,  the  multiplexing  factor,  the 
total  number  of  detectors,  and  the  number  of  pixels  can  be  calculated  as  a 
function  of  N„  for  a  given  flow  diameter,  total  number  of  fan  sources 
(fibers),  and  number  of  detectors  per  source.  The  results  for  several 
values  of  r,  Na,  and  n^  are  shown  in  Table  2.  As  shown  in  the  table  (and 
as  previously  noted),  w  is  nearly  independent  of  the  multiplexing  factor. 
Assuming  that  the  practical  upper  limit  on  the  diameter  and  atmospheric 
pressure  turbulent  flow  tube  is  4  in.,  noting  that  the  smallest  turbulence 
scale  size  for  this  flow  is  ~  1  mm  (Ref.  40),  requiring  Fm  to  be  an 
integer,  and  satisfying  all  the  constraints  mentioned  in  the  preceding 
discussion,  we  can  find  a  practical  geometry.  Such  a  geometry  is  shown  in 
Fig.  4.  In  this  figure,  Na  =  6,  0  =  30°,  and  F  =  3.86.  For  a  flow  diam¬ 
eter  of  4  in.,  96  fibers,  and  96  detectors  per  fan,  the  source-detecting 
ring  diameter  is  15.44  in.  (39  cm),  and  w  =  2.14  mm.  These  values  yield  a 
resolution  of  approximately  7000  pixels  for  a  multiplexing  factor  of  16  and 
a  total  of  576  detectors. 

B.  NONPLANAR  GEOMETRY 

In  the  previous  subsection,  we  considered  the  standard  fan  beam  geome¬ 
try,  where  the  source  and  detector  rings  had  equal  radii.  In  this  sub¬ 
section,  we  briefly  consider  a  nonplanar  geometry,  where  these  rings  are 
still  concentric  with  the  flow  but  have  different  diameters  (see  Figs.  5 
and  6).  There  are  three  disadvantages  to  a  nonplanar  design:  (1)  the 
reconstruction  of  the  field  from  the  projection  data  cannot  be  undertaken 
using  some  of  the  standard  algorithms,  (2)  a  right-angle  cross  section  of 
the  flow  is  no  longer  sampled,  and  (3)  the  instrument  is  optically  more 
complex.  These  drawbacks  may  be  outweighed  by  two  important  advantages. 
First,  the  detector  width  is  no  longer  dictated  solely  by  the  flow  diameter 
for  a  given  resolution.  Second,  the  radiation  reflected  from  the  trans¬ 
parent  flow  tube  from  the  fan  sources  is  not  directed  into  the  detectors. 
The  fraction  of  light  reflected  from  a  fan  is  considered  in  Section  III.B. 
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FLOW  DIAMETER (Inches)  =  4.00,  DETECTORS/FAN  =100.,  NUMBER  OF  FIBERS  =100. 
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Fig.  4 


Na  =  6 
e  =  30° 
F  =  3.86 


A  Practical  Geometry  for  Interrogating  a  Cylindrical  Cross  Sec¬ 
tion 
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Fig.  5.  Projected  View  of  a  Nonplanar  Fan  Beam  Geometry.  The  source  ring 
radius  Rg  is  smaller  than  the  detector  ring  radius  V 


Fig.  6.  A  Longitudinal  View  Along  a  Transparent  Flow  Tube  Depicting  the 
Displaced  Structure  of  the  Source  Ring  and  Detector  Ring  of  a 
Nonplanar  Fan  Beam  Instrument.  The  ray  emanating  from  the  upper 
ring  at  angle  <j>  is  partially  reflected  from  both  surfaces  of  the 
tube  before  emerging  from  the  opposite  side  of  the  tube  parallel 
to,  but  displaced  a  distance  5  from,  the  undeviated  path. 
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Consider  the  side  view  of  the  nonplanar  instrument  shown  in  Fig.  6.  A 
source  ring  of  radius  Rg  is  vertically  displaced  at  height  H  from  a  detec¬ 
tor  ring  of  radius  Rq.  It  is  straightforward  to  show  that  a  ray  of  light 
entering  the  side  of  the  tube  at  incidence  angle  <t>  is  displaced  upward  a 
distance 

6  =  T  sin«t>(  1  - )  (10) 

2  2 
/n  -  sin  4> 

where  T  is  the  thickness  of  the  tube  and  n  is  its  index  of  refraction.  For 
a  given  T,  Rs,  RD,  and  r,  4>  can  be  found  as  a  function  of  H  by  requiring 
the  ray  to  strike  the  detector  ring,  i.e. 

H  =  A1  -  6  +  A2  -  6  +  A3  (11) 


where 


A.  =  (R  -  r)tan<|> 

I  S 


A2  =  2r  tan<|> 


A^  =  (Rp  -  r)tan<(> 


or 


(R  +  Rn  -  2r)tan$  -  2Tsin<|>  (l  - 

S  u 


cost > 


2  2 
/n  -  sm  i 


:)  -  H  =  0 


(12) 


In  general,  one  wishes  to  choose  H  as  small  as  possible,  consistent  with 
the  requirement  that  the  reflected  rays  from  the  inner  or  outer  surfaces  of 
the  tube  do  not  strike  the  detector  ring.  "This  requires  x1  *  RD  and 

x2  *  RD»  where 


x,  = 


f1 


1  tan<t> 


+  r  +  T 


x„  = 


2  "  tan i)> 


+  r  +  T 


(13) 

(14) 
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y1  =  H  -  (Rg  -  r  -  T)tan4> 


(15) 


y2  =  y1  -  2T  tan<j>'  ( 16) 

4>  =  sin-1  (17) 

Moreover,  we  require  that  the  ray  miss  the  detector  by  more  than  half  the 
fan  thickness,  tf.,  plus  half  the  detector  height,  hD,  or 

t1  =  (x1  -  RD)tan$  >  (tf  +  hD)/2  (18) 

t2  =  (x2  -  RD)tan4>  >  (tf  +  hD)/2  (19) 


The  procedure  is  to  choose  H  as  small  as  possible,  keeping  it  greater  than 
tf.  +  hD  so  the  beam  does  not  strike  the  opposite  side  of  the  source  rings; 
to  solve  Eq.  (12)  for  4>  for  the  chosen  H;  then  to  check  that  inequalities 
in  Eqs.  (18)  and  (19)  are  satisfied. 

Figure  7a  displays  the  solution  of  Eq.  (12)  for  r  =  2.0  in.,  T  = 

0.125  in.,  Rd/Rs  =  3,  and  n  =  1.5  for  Rg/r  =  2, 3, 4, 6, 8.  Figure  7b  plots 

ti»  t2,  and  (tf  +  hD)/2  for  these  values,  letting  tf  =  1  cm  and  hs  = 

1  cm.  Figure  7b  indicates  that  the  reflected  rays  intersect  the  detector 

ring  between  5  and  6°.  Thus,  choosing  a  geometry  where  <t>  is  less  than  4° 

ensures  that  reflected  radiation  will  not  enter  the  detectors.  These  plots 

are  not  functions  of  R_/r. 

s 
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NONPLANAR  GEOMETRY  R  /r  =  g 


2 

4  6  8 

TILT  ANGLE  (DEGREES) 

REFLECTED  RAY  QUANTITIES 


(t,  +  hA)/2 


TILT  ANGLE  (DEGREES) 


(a)  The  Relationship  Between  the  Tilt  Angle  <(>  shown  in  Fig.  6  and 
the  Ring  Separation  in  the  Nonplanar  Design;  (b)  Plot  of  the 
Reflection  Quantities  Discussed  in  Text.  Only  when  the  dashed 
curve  intersects  the  solid  curve  does  the  reflected  radiation 
from  the  first  two  surfaces  of  the  tubes  strike  the  detectors. 


III.  TUBE  PERTURBATIONS 


A.  OVERVIEW 

One  of  the  more  important  applications  of  the  proposed  tomographic  in¬ 
strument  is  the  interrogation  of  a  flow  confined  in  a  cylindrical  tube. 

The  observation  port  should  provide  360°  access  to  the  cross  section  of  the 
flow  at  the  position  of  interest.  Ideally,  the  view  port  will  be  a  short, 
transparent  tube  with  the  optical  quality  of  an  ophthalmic  lens.  This  view 
port  will  be  a  glass,  quartz,  or  sapphire  thin-walled  cylinder  with  an 
inner  diameter  equal  to  the  diameter  of  the  flow  (~  5-15  cm),  a  length  of  1 
to  2  cm,  and  a  wall  thickness  of  2  to  6  mm.  The  view  port  will  be  posi¬ 
tioned  within  the  break  of  the  insulated  flow  tube.  The  length  of  the  port 
will  be  kept  as  small  as  possible  to  minimize  thermal  losses  from  the  gas. 

Figure  8  depicts  a  fan  of  radiation  striking  such  a  view  port.  A  por¬ 
tion  of  each  ray  within  the  fan  is  reflected  from  the  outer  surface  of  the 
tube.  In  a  planar  instrument,  this  reflected  radiation  strikes  one  of  the 
detectors  on  the  source-detection  ring.  The  portion  of  the  ray  that  passes 
through  the  first  surface  continues  to  be  divided  and  reflected  around  the 
tube,  as  shown  in  Fig.  8b  for  the  ray  launched  at  angle  a.  The  transmitted 
intensities  marked  L^L^Ly,  etc.,  escape  the  tube  and,  for  the  planar 
geometry,  also  strike  the  detectors.  The  transmitted  intensity  marked 
is  the  probe  beam  (the  heavy  line),  which  suffers  a  deviation  as  it  passes 
through  the  four  interfaces.  The  arc  length  displacement  of  this  beam  on 
the  detector  ring  from  its  unperturbed  position  is  marked  As.  The  unper¬ 
turbed  trajectory  (the  dashed  line)  follows  a  different  path  through  the 
gas  than  does  the  light  ray.  The  mathematical  algorithms  used  to 
reconstruct  a  cross  section  of  the  flow  assume  an  undeviated  path.  The 
intensity  transmitted  through  the  first  interface  t-j  may  be  totally 
reflected  at  the  second  interface  if  the  incidence  angle  at  this  interface 
is  greater  than  the  critical  angle  (=  41.8°  for  an  index  of  refraction 
ratio  of  1. 5/1.0).  The  deviation  of  the  ray  is  discussed  in  Section 
II I. B.  The  problem  of  the  reflected  and  transmitted  intensities  is 
considered  in  Section  III.C. 
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Fig.  8.  (a)  Several  Rays  of  a  Fan  Striking  the  Outside  of  a  Transparent 

Flow  Tube;  (b)  The  Fate  of  a  Single  Ray  at  Angle  a  in  a  Fan  as  It 
Traverses  the  Tube  and  Strikes  the  Detector  Ring. 
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B.  RAY  DEVIATION 


The  view  port  is  equivalent  to  two  cylindrical  lenses  separated  by  a 
distance  equal  to  the  sum  of  their  outer  (R^)  and  inner  (R2)  radii  of  cur¬ 
vature.  Paraxial  ray  theory  can  be  used  to  calculate  the  focal  length  of 
the  left  side  of  the  tube 


f  =  (n  -  D(R^  -  (20) 

where  R^  (R2)  is  the  radius  of  curvature  of  the  first  (second)  surface 
struck  by  the  rays.  The  radius  is  positive  if  the  center  of  curvature  lies 
to  the  right  of  the  surface  and  negative  if  it  lies  to  the  left  of  the 
surface.  Accordingly,  the  radii  are  positive  for  the  left  side  of  the  tube 
and  negative  for  the  right  side  of  the  tube.  The  net  effect  is  that  both 
sides  of  the  tube  act  as  negative  (diverging)  lenses.  For  example,  if  | R ^  | 
=  5.5  cm  and  |R2|  =  5.0  cm,  f1  =  fg  =  -110  cm.  For  two  identical  thin 
lenses  separated  by  a  distance  d,  the  front  and  back  focal  lengths  are 
identical  and  are  given  by  (Ref.  Ml) 


f (d  -  f) 
fb  "  d  -  2f 


(21) 


If  d  =  R1  +  Rg,  ffb  =  -57.5  cm.  Thin  lens  theory  cannot  be  used  to  calcu¬ 
late  correct  deviations,  even  for  thin-wall  cylinders,  because  the  geom¬ 
etries  under  consideration  require  large  fan  angles  (nonparaxial  rays).  An 
exact  ray  trace  calculation  is  required. 

Consider  the  geometry  shown  in  Figs.  9a  and  10.  A  ray  a  distance  di 
from  the  outer  surface  of  the  tube  is  launched  at  an  angle  a.,  which  is  a 
distance  y.^  above  the  horizontal.  This  ray  strikes  the  tube  surface  at  a 
point  a  height  y^  above  the  horizontal,  making  an  angle  We  wish  to 
calculate  y1  and  B1  as  a  function  of  R^ ,  cu,  yi?  and  di.  We  note 


y  1  =  (c  +  d +  s.|  )tan 


(22) 


31 


32 


-  R^  cos  6^ 

c  =  y . /tan  a . 

J 1  1 

y 1  =  R1  sin  e1 

Substituting  Eqs.  (23) -( 25 )  into  Eq.  (22)  and  solving,  we  find 


(23) 

(24) 

(25) 


where 


. -1  ♦  /t!  ■  t!  t  I, 

S,  =  cos  ( - j - - - ) 

t  +  1 


k  =  [y±  +  (d.^  +  R1)t]/R1 


(26) 


(27) 


t  =  tan  cu 

and  y1  is  given  by  Eq.  (25). 

Referring  to  Fig.  9b  for  the  second  surface,  we  wish  to  find  b2  and 
y2,  given  y1 ,  ,  B^ ,  R.,,  and  R2.  Here 


y2  =  (a  -  s1  +  R^  -  R2  +  s2)tan 

(28) 

s2  =  R2(1  -  cose2) 

(29) 

a  =  y^/tana^ 

(30) 

y2  =  R2  sin  b2 

(3D 

Substituting  Eq.  (23)  and  Eqs.  (29M31)  into  Eq.  (28),  we  find  &2  is  given 
by  Eq.  (26),  where 


k  =  (y 1  +  R1  t  cosb1)/r2 

(32) 

t  =  tana^ 

and  y2  is  given  by  Eq.  (31). 

Surface  3  is  shown  in  detail  in  Fig.  9c,  where  the  unknowns  B^  210(1  y3 


are  obtained  in  terms  of  a 2,  y2,  R2,  and  B2 

y^  =  [(a  -  s2)  +  2R2  -  s^ltanag  (33) 

a  =  y2/tanc*2  (34) 

S3  =  R2(1  -  cosB3)  (35) 

y3  =  R2  sinB^  (36) 


Substituting  Eq.  (29)  and  Eqs.  (34)-(36)  into  Eq.  (33),  we  find  B^  is  given 
by  Eq.  (26),  where 


k  =  yD/R_  -  t  cosBo 

*  (37) 

t  =  -tana2 

and  y3  is  given  by  Eq.  (36).  The  geometry  for  surface  4  is  shown  in  Fig. 

9d 


y4  =  (a  +  s^  +  R1  -  R2  -  s^tanct^ 
a  =  y^/tano^ 
s4  =  R-,0  -  cosB4) 
y4  =  R1  sinB4 


(38) 

(39) 

(40) 

(41) 
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Substituting  Eq.  (35)  and  Eqs.  ( 39 ) - (4 1 )  into  Eq.  (38),  we  find  6^  is  given 
by  Eq.  (26),  where 


k  =  (y^  +  R2t  cosb^/R.] 


(42) 


t  =  -tano^ 

and  yjj  is  given  by  Eq.  (41). 

The  composite  drawing  is  shown  in  Fig.  10.  Each  angle  o  must  be 
related  to  the  preceding  a  through  Snell's  law 


n.  sine.  =  n,  sine.  (43) 

1  1C  w 

where  e.  (e.  )  is  an  incidence  (transmitted)  angle  measured  with  respect  to 

1  t 

the  surface  normal  and  n  is  the  index  of  the  medium.  Then  for  surface  1 


6i1 


ai  +  B1 


-1  n1 

).  ,  =  sin-  ( —  sin  e. , ) 
tl  ln2  i1' 


6t1  -  61 


(44) 


for  surface  2 


(45) 


for  surface  3 
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(46) 


and  for  surface  4 


(47) 


Finally,  we  note  that  the  ray  strikes  the  detector  ring  at  a  height  yf 


yf  =  Rd  sinBf 


(48) 


where  Bf  is  given  by  Eq.  (26)  and 

k  =  (y^  +  R-] t  cosB^)/Rd 


(49) 


t  —  “  t  Bn  ct|| 

i 

The  ray  would  have  struck  the  detector  ring  at  height  y~  if  the  tube  had 

i  * 

not  been  present,  where  yf  is  given  by 


i  i 

yf  =  Rd  sinBf 


(50) 


and  is  given  by  Eq.  (26),  with 


k  =  [yt  -(d.  +  R -j ) t ] /RD 
t  =  -tan  cu 

Finally,  the  arc  length  difference  is 


(51) 
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As  =  RD(e^  -  sp 


(52) 


Figure  11  shows  the  are  length  deviation  [Eq.  (52)]  as  a  function  of 
angle  within  the  fan  for  Na  =  6  for  three  tube  inner  diameters  and  three 
wall  thicknesses.  Table  3  lists  the  undeviated  height  [Eq.  (50)],  the  ray- 
refracted  height  [Eq.  (48)],  and  the  arc  length  difference  [Eq.  (52)], 
where  Na  =  6,  D  =  4  in.,  and  t  =  1/8  in.  for  50  angles  within  a  fan.  These 
results  clearly  indicate  that  the  deviations  are  large  compared  to  the 
detector  widths  for  all  but  the  most  central  angles.  Also,  it  is  desirable 
to  use  as  thin  a  wall  as  possible,  and  for  a  given  tube  thickness,  as  large 
an  inner  tube  diameter  as  possible.  This  is  because,  for  a  given  value  of 
Na,  nd,  and  nj.,  the  space  available  for  the  detector  increases  linearly 
with  the  flow  diameter,  whereas  the  deviations  increase  less  than  linearly 
with  R2.  The  expressions  derived  above  can  be  used  to  preprocess  the  pro¬ 
jection  data  before  reconstruction  of  the  cross  section  of  flow  or  can  be 
used  to  modify  the  reconstruction  algorithm  itself.  Given  the  magnitude  of 
those  deviations,  it  may  be  necessary  to  limit  the  fan  beam  to  the  central 
three  quarters  of  the  flow  and  to  use  an  onion  peeling  algorithm  to  help 
reconstruct  the  flow. 

As  noted  in  Section  III. A,  because  the  actual  ray  path  and  the  undev¬ 
iated  ray  path  (dashed  line  in  Fig.  10)  are  not  the  same  inside  the  tube, 
an  incorrect  reconstructed  field  may  result.  The  reconstruction  algorithms 
assume  that  the  undeviated  path  is  taken.  Correcting  the  algorithm  only 
for  the  displacement  As  at  the  detector  ring  is  acceptable  only  if  the 
difference  between  the  undeviated  path  and  the  actual  path  is  less  than  or 
equal  to  the  resolution  (detector  spacing). 

Consider  the  difference  between  the  coordinates  of  the  end  points  of 
the  two  rays  at  the  points  where  they  intersect  the  inner  wall  of  the 
tube.  The  coordinates  of  intersection  of  the  undeviated  ray  with  the  inner 
circle  of  radius  R2  centered  at  the  origin  are 
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LAUNCH  ANGLE  (DEGREES) 


Fig.  11.  Deviation  at  Detector  Ring  of  Ray  Passing  Through  a  Tube  for 
Three  Different  Inner  Diameters  and  Three  Wall  Thicknesses 
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Table  3.  Refraction  Calculations  for  Tubes 
D  =  4.0  in.,  t  =  1/8  in.,  Na  =  6 


unch  anqlef  rad) 

undeviated  ht(cra) 

refracted  ht(cra) 

deviation(mm) 

0.0052408 

0.20547 

0.2215? 

0.16102 

0.0104816 

0.41110 

0.44354 

0.32452 

0.0157224 

0.61654 

0.66526 

0.48746 

0.0209632 

0.82198 

0.88712 

0.65205 

0.0262040 

1.0273 

1.1089 

0.81763 

0.0314448 

1.2325 

1.3308 

0.98560 

0.0366856 

1.4375 

1.5528 

1.1560 

0.0419264 

1.6424 

1.7748 

1.3293 

0.0471672 

1.8471 

1.9970 

1.5058 

0.0524080 

2.0517  " 

2.2193 

1.6860 

0.0576488 

2.2560 

2.4417 

1.8707 

0.0628896 

2.4600 

2.6643 

2.0605 

0.0681304 

2.6638 

2.8871 

2.2553 

0.0733712 

2.8673 

3.1101 

2.4566 

0.0786120 

3.0705 

3.3334 

2.6646 

0.0838528 

3.2733 

3.5569 

2.8802 

0.0890936 

3.4758 

3.7808 

3.1040 

0.0943344 

3.6779 

4.0051 

3.3373 

0.0995752 

3.8796 

4.2299 

3.5808 

0.1048160 

4.0808 

4.4552 

3.8355 

0.1100568 

4.2817 

4.6811 

4.1031 

0.1152976 

4.4820 

4.9077 

4.3844 

0.1205384 

4.6819 

5.1351 

4.6811 

0.1257792 

4.8812 

5.3634 

4.9952 

0.1310200 

5.0800 

5.5927 

5.3284 

0.1362608 

5.2782 

5.8233 

5.6828 

0.1415016 

5.4759 

6.0552 

6.0611 

0.1467424 

5.6730 

6.2888 

6.4662 

0.1519832 

5.8694 

6.5242 

6.9015 

0.1572240 

6.0652 

6.7617 

7.3707 

0.1624648 

6.2603 

7.0017 

7.8788 

0.1677056 

6.4548 

7.2447 

8.4308 

0.1729464 

6.6485 

7.4910 

9.0334 

0.1781872 

6.8415 

7.7413 

9.6944 

0.1834280 

7.0338 

7.9963 

10.423 

0.1886688 

7.2252 

8.2569 

11.231 

0.1939096 

7.4159 

8.5242 

12.133 

0.1991504 

7.6058 

8.7996 

13.147 

0.2043912 

7.7948 

9.0848 

14.296 

0.2096320 

7.9830 

9.3821 

15.611 

0.2148728 

8.1703 

9.6946 

17.132 

0.2201136 

8.3567 

10.026 

18.914 

0.2253544 

8.5422 

10.383 

21.037 

0.2305952 

8.7267 

10.773 

23.615 

0.2358360 

8.9103 

11.209 

26.827 

0.2410768 

9.0929 

11.712 

30.971 

0.2463176 

9.2745 

12.319 

36.588 

0.2515585 

9.4551 

13.102 

44.829 

0.2567993 

9.6347 

14.260 

58.951 
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2 

-m.  a  + 
1 


(53) 


x 


2,3 


2 

m.  +  1 

l 


l2)  +  R 


2 

2 


where  x^,  yj,  (x^,  y^)  are  the  coordinates  of  the  left  (right)  points  of 
intersection  in  Fig.  10,  m^  =  tancu,  and  1  =  c  +  d^  +  R^.  The  y  coor¬ 
dinates  of  the  intersection  of  the  actual  ray  at  the  inner  tube  radius  are 
given  by  Eqs.  (31)  and  (36).  The  x  coordinates  are  then 


2,3 


(y2,3  "  b1,2)/m1,2 


(54) 


where  b1>2  =  y-|>2  "  m1,2x1,2>  m1,2  =  tan“i  2’  X1  =  ^y1  “  bi^/mi*  bi  = 
y2  -  nrxi,,  y1  is  given  by  Eq.  (25),  and  2  are  given  in  Eqs.  (44)  and 

(45). 

Figure  12  displays  the  results  of  the  calculations  for  a  4.0  in.  diam¬ 
eter  flow  and  a  tube  thickness  of  0.125  in.  The  quantities  (x„  ~  -  x'  _), 
r  2  2i1/2  ^ >  5  c- , 5 

(y2>3  "  y2  3^’  Ux2  3  ~  x2  3^  +  ^y2  3  “  y2  3^  J  *  and  the  normalized 
position  of  the  crossing  point  of  the  two  rays  (xc  -  x2)/(x^  -  x2),  (yc  - 

y2)/(y3  -  y2)  are  shown  in  Figs.  12a,  12b,  12c,  and  12d,  respectively.  The 
crossing  coordinates  are  given  by 


x  = 


(b.  -  b2)/(m2 


mi} 


y„  = 


m„  x  +  b„ 
2  c  2 


(55) 


These  results  show  that  the  deviations  at  the  end  points  approach  twice  the 
detector  spacing.  This  effect  could  degrade  the  resolution  by  a  factor 
of  ~  2  near  the  inner  wall  of  the  tube  if  not  corrected  in  the  reconstruc¬ 
tion  algorithm. 
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Fig.  12.  Curves  Depicting  the  Difference  Between  the  Undeviated  Ray  and 
the  Actual  Ray  Inside  the  Tube,  (a)  x  coordinate  differences  at 
the  ray  intersection  points  with  the  tube  inner  wall;  (b)  y  co¬ 
ordinate  differences;  (c)  distance  between  the  end  points  of  the 
two  rays  at  the  inner  wall  intersection;  (d)  normalized  crossing 
point.  A  value  of  0.5  indicates  crossing  at  the  midpoint  between 
the  intersection  of  the  actual  ray  with  the  inner  wall. 
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REFLECTION  LOSSES 


The  intensities  of  the  reflected  and  transmitted  rays  shown  in  Fig.  8 
can  be  calculated  using  Fresnel's  equations  (Ref.  41) 


sin(6.  -  e  ) 

R ,  ( e . ,  e.  )  =  - rM 

1  i’  t  ‘sin  0.  +  0. )j 

l  t 


(56) 


tan(e.  -  et) 


R||  (V  et)  *  ItSje—  e j) 

1  t 


(57) 


23100,.  cos0.  n,.  cose. 


WWW  •  11,  WWW  w  ■ 

T,  (0.,  e. ,  n.,  n. )  =  [ T-- - — ±1  [— - -) 

l  i’  t’  i’  t  lsin(e.  +  0  ) J  '•n.  cos0.; 

1  w  XI 


(58) 


2  sine.  cose. 

C  1 


nt  cos0t- 


|  |  ^  8i ’  0t’  ni’  nt^  "  ^sin(e.  +  0.  )  cos(6.  -  e.)l  ^n.  cos0.^ 

1  t  Xu  XI 


(59) 


where  R^  (R||)  and  (T 1 1 )  are  the  reflected  and  transmitted  intensities, 
respectively,  for  s-polarized  (p-polarized)  light.  The  incident  ray  in  a 
medium  with  index  n.^  has  an  incident  angle  e^  is  reflected  from  the  sur¬ 
face  with  angle  0^,  and  is  transmitted  into  the  medium  with  index  nfc  at 
angle  6^ •  All  angles  are  measured  with  respect  to  the  surface  normal.  The 
angles  are  related  through  Snell's  law  [Eq.  (43)]  and  the  law  of  reflec¬ 
tion,  e.  =  e  . 

’  l  r 

Consider  the  geometry  shown  in  Fig.  13.  From  A ABC,  6^  =  0^  -  e  ; 
ABCD,  fig  =  6^2  “  0i3»  ACDE,  6^  =  “  0i3’  ACEF,  5^  =  0.^  -  0^,.;  etc. 

From  AABD,  0^  =  0.^;  ADEF,  e^  =  0^;  etc.  Thus,  we  conclude 


r  1 


=  e 


13 


=  e.c 
i5 


=  0i7 


=  0. 


10 


(60) 
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Fig.  13.  Geometry  and  Definition  of  Parameters  Used  to  Calculate  the 
Intensities  of  the  Reflected  Radiation 


(61) 


ei2  =  ei4  =  0i6  =  0i8  =  1  6ie 


61  =  S2  ~  S3  =  64  =  -  =  6  =  e2  -  e1 


Because  0^,  02,  and  0^  are  known  from  the  previous  subsection  [Eqs. 
(26),  (27),  (32),  and  (45)],  0ie  and  6  are  known.  Thus,  only  0iQ  must  be 
found  to  have  all  angles  of  incidence.  Angle  0^q  is  obtained  from  0^. 
Applying  the  cosine  law  twice  to  ABCD,  we  have 

—  2  2 

BD  =  +  R2  —  2R^2  cos6 


R2  =  (BD)2  +  R2  -  2  (BD)  R1  cos0.3 


Solving  for  0i3,  we  find 


e.  =  e._  =  cos 
10  13 


-  R2  cosS 

2  2  1/2 
(R*  +  R*  -  2R1R2  cosS) 


Finally,  the  angles  of  transmission  (refraction)  are  obtained  from  Eq.  (43) 


.  -1  2  . 

0.  =  sin  —  sin  0. 

te  ie 


-1  2 

0.  =  sin  —  sin  0. 

to  n^  10 


If  the  intensity  of  the  incident  ray  is  assigned  a  normalized  value  of 
1,  then  the  intensities  of  the  beams  labeled  Ln,  where  n  =  1,2,3,  etc.,  are 


L1  =  R1 


L2  =  T1Te 


S  =  T1ReTo 


L2n  =  {Wn~lL2 


L0  ,  =  (R  fi  )n_1L- 
2n+1  e  o  3 


(65) 


In  these  relations,  n  refers  to  the  surface  encountered,  and  we  used 


T.  =  Tc  =  T,  =  . . . .  =  T 

3  o  7  o 


(66) 


The  signal  beam  is  intensity  L2.  All  other  values  of  Ln  are  losses 
from  this  signal  beam.  The  total  exterior  and  interior  losses  are 


n=0  L2lw1 


(67) 

(68) 
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Finally,  the  intensity  of  a  beam  traversing  the  entire  tube  (the  projection 
of  L2  through  the  second  wall)  is 


L' ' 
L2 


=  L2T3T4 


(69) 


where  and  are  Eqs.  (58)  and  (59)  evaluated  for  the  incident  and 
transmitted  angles  at  the  back  two  interfaces  of  the  tube  wall.  Total 
internal  reflection  occurs  at  dense-to-rare  interfaces  when  the  angle  of 
the  transmitted  beam  is  >  90°,  that  is,  when 


or 


n 

sin  6.  >  — 

le  -  n2 


n 


1 


sine.  >  — 
10  -  n2 


(70a) 

(70b) 


Fan  beam  geometries  that  satisfy  either  Eqs.  (70a)  or  (70b)  must  be 
avoided . 

Equations  (56 )-(69)  can  be  used  to  calculate  all  relevant  intensities 
in  the  fan  beam.  Figures  14-18  exhibit  the  results  of  some  of  these  cal¬ 
culations.  Figures  14-15  display  the  intensities  of  the  transmitted 
beam  (LJJ),  the  interior  transmitted  beam  (L2),  the  first  surface  reflection 
loss  (L.|),  the  first  transmitted  loss  out  of  the  tube  (L^),  the  total 
external  loss  [Eq.  (67)],  and  the  total  internal  loss  [Eq.  (68)],  for  s-, 
p-,  and  unpolarized  beams  as  a  function  of  ray  angle  within  the  fan  beam 
for  N  =  6;  an  inner  diameter  of  4  in.;  and  tube  wall  thicknesses  of  1/8  and 
1/4  in.  Figures  16  and  17  exhibit  the  same  calculations  for  a  tube  with  an 
inner  diameter  of  6  in.  and  wall  thicknesses  of  1/8  and  1/4  in.  Figure  18 
displays  results  for  a  6  in.  tube  with  N  =  3. 

From  these  curves,  we  can  draw  a  number  of  conclusions.  The  first 
surface  reflection  loss  is  greater  for  the  6  in.  tube  than  for  the  4  in. 
tube  (13£  vs  1 0/6 )  and  is  greater  for  the  1/8  in.  wall  than  for  the  1/4  in. 
tube  ( 10%  vs  1 3%  at  the  maximum  ray  angle).  The  transmitted  beam  intensity 
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PERCENT  PERCENT  PERCENT 


LAUNCH  ANGLE  (DEGREES) 


Fig.  14.  Reflected  and  Transmitted  Intensities  for  All  Ray  Angles  Within  a 
Fan  for  Flow  Diameter  of  4  in.,  Tube  Wall  Thickness  of  0.125  in., 
and  N„  =  6 

a 
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Fig.  15.  Reflected  and  Transmitted  Intensities  for  all  Ray  Angles  Within  a 
Fan  for  Flow  Diameter  of  4  in..  Tube  Wall  Thickness  of  0.25  in., 
and  Na  =  6 
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16.  Reflected  and  Transmitted  Intensities  for  All  Ray  Angles  Within  a 
Fan  for  Flow  Diameter  of  6  in.,  Tube  Wall  Thickness  of  0.125  in., 
and  N.  =  6 

a 
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Fig.  17.  Reflected  and  Transmitted  Intensities  for  All  Ray  Angles  Within  a 
Fan  for  Flow  Diameter  of  6  in..  Tube  Wall  Thickness  of  0.25  in., 
and  Na  =  6 
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Fig.  18.  Reflected  and  Transmitted  Intensities  for  All  Ray  Angles  Within  a 
Fan  for  Flow  Diameter  of  6  in..  Tube  Wall  Thickness  of  0.25  in., 
and  Na  =  3 
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is  slightly  better  for  the  4  in.  tube  than  for  the  6  in.  tube  (by  1  to  2% 
at  high  ray  angles)  and  is  slightly  better  for  the  1/4  in.  wall  than  for 
the  1/8  in.  wall.  The  first  transmitted  loss  out  of  the  tube  (Lg)  is  some¬ 
what  greater  for  the  4  in.  tube  and  is  greater  for  the  1/4  in.  wall  (35%  vs 
28%) .  The  total  external  loss  is  the  same  for  both  tube  diameters  and  very 
slightly  greater  for  the  1/8  in.  tube  than  for  the  1/4  in.  tube.  The  total 
internal  loss  is  greater  for  the  6  in.  diameter  tube  than  for  the  4  in. 
diameter  tube  (3%  vs  ~  2%)  and  is  greater  for  the  1/8  in.  wall  than  for  the 
1/4  in.  wall  (4J  vs  3%).  No  significant  difference  is  seen  between  the 
N  =  3  and  N  =  6  calculations.  In  general,  the  4  in.  diameter  tube  is 
slightly  better  than  the  6  in.  diameter  tube,  and  the  1/4  in.  wall  is 
slightly  better  than  the  1/8  in.  wall. 

These  curves  indicate  that  it  would  be  significantly  better  to  use 
p-polarized  light.  Unfortunately,  polarization-preserving  fibers  currently 
exist  only  for  single-mode  optic  fibers.  Because  of  the  high  intensities 
experienced  during  a  short  (nanosecond)  laser  pulse,  only  multimode,  macro¬ 
fibers  are  relevant  to  this  work. 

In  summary,  the  calculations  from  this  subsection  and  the  previous 
subsection  favor  a  thin-walled,  4  in.  diameter  flow  tube  view  port. 
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IV.  FAN  BEAM  DISTORTION  CAUSED  BY  THE  GAS 


Variation  of  the  index  of  refraction  of  the  gas  causes  distortion  in 
the  fan  beam.  The  index  is  a  function  of  the  gas  density  which,  in  turn, 
depends  on  the  temperature  and  pressure.  For  flow  tube  applications,  pres¬ 
sure  gradients  are  expected  to  be  small,  and  temperature  variations  cause 
the  principal  distortions  of  the  rays.  In  this  section,  we  consider  the 
deviation  of  a  ray  in  the  fan  beam  due  to  a  cylindrically  symmetric  temper¬ 
ature  variation. 

The  path  of  a  ray  traversing  a  medium  with  a  continuously  varying 
index  of  refraction  is  given  by  (Ref.  42) 


dr 

ds 


)  =  vn 


(71) 


where  r  is  the  position  vector  along  the  ray  path  and  s  is  the  path 
length.  If  we  refer  to  the  geometry  shown  in  Fig.  19a,  Eq.  (71)  can  be 
written 


e<yi>  =  jj-  f 


o  -00 


36 

3y 


dx 


(72) 


where  n  =  1  +  6,  nQ  is  the  index  of  the  unperturbed  medium,  —  =  0(yi),  and 

it  is  assumed  that  the  variation  of  the  index  is  small  compared  to  1 

(paraxial  approximation).  This  approximation  is  always  true  for  gases  at 

9  5  96 

or  below  atmospheric  pressure  and  allows  ds  =  dx  and  —  =  — .  If  we  trans¬ 
form  Eq.  (72)  to  polar  coordinates,  assuming  6  is  a  function  of  r  only 

(cylindrical  symmetry),  noting  that  the  integrand  is  an  even  function  of  x 


3  6 

and  that  —  =  0  for  r  >  R,  we  have 
9r 


e(yi) 


n 


o 


; 

*i 


Mill 

3r 


(73) 
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Fig.  19.  Definition  of  Quantities  Used  to  Calculate  the  Ray  Deviation 
Caused  by  a  Cylindrically  Symmetric,  Continuously  Varying  Index 
of  Refraction 
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For  a  fan  beam  geometry,  we  wish  to  find  e(or  ),  where  or  is  the  angle 
of  ray  (i)  in  the  fan  (Fig.  19b).  The  relationship  between  yi  and  or  is 
especially  simple  for  cylindrical  symmetry,  because  each  value  of  or  cor¬ 
responds  to  a  new  yA.  Figure  19c  indicates  that  a  ray  launched  at 
angle  or  in  Fig.  19b  has  exactly  the  same  geometry  as  Fig.  19a  with  a  new 
coordinate  system  (x * ,y ' )  if  the  flow  has  cylindrical  symmetry.  According¬ 
ly 


yi  =  xftancr  (74) 

The  index  of  refraction  of  a  gas  is  a  function  of  the  density,  i.e. 

(n-1)=6«p 


For  constant  pressure  and  an  ideal  gas 


where  C  is  a  constant  and  T  is  the  absolute  temperature  in  Kelvin.  Noting 
fiNPT  =  0.000293,  we  have  C  =  8.79  x  10”2  K.  If  the  temperature  at  the 
center  of  the  flow  is  Tm,  the  edge  temperature  is  TQ  and  we  allow  the  tem¬ 
perature  to  vary  as  a  power  of  r  between  these  values,  i.e. 


T(r )  =  Tffl  -  krn  n  =  1,2,3,... 

where  ( 75 ) 

k  =  (T  -  T  )/Rn 
m  0 


Then 


6(r)  =  C(T  -  kr11)-1 
in 

and  the  gradient  of  the  index  becomes 


(76) 
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(77) 


2 

d6(r)  .  (n—  1 )  .  n, 

- v  —  =  nkCr  /(T  -  kr  ) 

dr  m 

Using  Eqs.  (77)  and  (74)  in  Eq.  (73),  and  numerically  integrating  Eq. 
(73),  we  can  calculate  the  deviation  as  a  function  of  c^,  geometry,  TQ,  Tro, 
and  temperature  gradient.  The  results  of  these  calculations  for  several 
geometries  and  for  five  powers  of  r  (linear  through  r^)  are  given  in  Figs. 
20  and  21.  These  results  indicate  that,  in  all  cases,  the  deviations  at 
the  detector  ring  caused  by  the  temperature  gradients  are  small  compared  to 
the  detector  width  (w  =  2  mm  for  D  =  4  in.  and  w  =  3  mm  for  D  =  6  in.). 

The  curves  also  indicate  that  the  deviations  increase  with  increasing  N„ 

a 

and  D.  This  increase  is  expected  because  the  detector-ring  diameter 
increases  with  both  of  these  parameters.  The  deviations  are  very  sensitive 
to  the  difference  between  the  center  and  edge  temperatures.  In  Fig.  21, 
changing  the  edge  temperature  from  300  to  500  K  reduces  the  maximum  devia¬ 
tion  by  nearly  a  factor  of  3.  In  all  curves,  in  the  angular  central  two- 
thirds  of  the  tube,  the  deviation  is  nearly  independent  of  the  functional 
form  of  the  temperature  gradient.  Thus,  if  the  peak  and  edge  average  tem¬ 
peratures  can  be  estimated,  average  corrections  can  be  applied  to  this 
region,  either  a  priori  or  iteratively. 

This  calculation  can  also  be  used  to  estimate  the  ray  deviation  caused 
by  a  turbulence  cell.  Figure  22  exhibits  the  deviation  at  the  detector 
ring  due  to  a  2  mm  diameter  bubble  with  a  center-to-edge  density  variation 
of  a  factor  of  2.  At  high  angles  and  for  radial  variations  of  ~  r2-r^, 
the  deviation  is  approximately  0.05  mm.  This  value  is  not  a  strong  func¬ 
tion  of  the  size  of  the  structure.  A  nonrandom  cumulative  effect  of  sev¬ 
eral  such  structures  could  be  appreciable  and  might  restrict  the  degree  of 
turbulence  tolerable  for  high-resolution  imaging. 

Last,  we  note  that  the  index  of  refraction  of  a  medium  changes  rapidly 
near  an  absorption,  as  shown  in  Fig.  23.  The  refractive  index  n  varies 
from  less  than  1  to  greater  than  1  as  radiation  is  scanned  across  an 
absorption  curve.  This  is  not  a  problem  for  frequency-stable  radiation 
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Fig.  20.  Deviation  at  the  Detector  Ring  of  a  Ray  Passing  Through  a  4  in 
Diameter  Flow  with  a  Cylindrically  Symmetric  Variation  in  Tem 
perature  for  Five  Different  Temperature  Gradients 
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Fig.  21.  Deviation  at  the  Detector  Ring  of  a  Ray  Passing  Through  a  6  in. 

Diameter  Flow  with  a  Cylindrically  Symmetric  Variation  in  Tem¬ 
perature  for  Five  Different  Temperature  Gradients 
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Fig.  22.  Deviation  at  the  Detector  Ring  of  a  Ray  Passing  Through  a  2  mm 
Turbulence  Bubble  with  a  Cylindrically  Symmetric  Variation  in 
Density  for  Five  Different  Density  Gradients.  The  edge  density 
is  a  factor  of  2  greater  than  the  center  density. 


(c) 


Fig.  23.  Ideal  Dispersion  Curves  for  an  Absorption  with  Different  Amounts 
of  Damping  and  Absorption.  (a)  Strong  absorption,  strong  damp¬ 
ing;  (b)  strong  absorption,  weak  damping;  (c)  weak  absorption, 
strong  damping;  (d)  weak  absorption,  weak  damping.  The  absorp¬ 
tion  coefficient  is  K. 


with  a  bandwidth  much  less  than  the  absorption  linewidth  centered  on  the 
peak  absorption.  Any  drift  in  wavelength,  however,  could  cause  appreciable 
refraction  if  the  probe  gas  is  a  major  species.  This  is  not  expected  to  be 
a  problem  if  the  concentration  of  the  species  is  less  than  0.1? 


V.  DATA  REDUCTION:  LINEARITY,  DYNAMIC  RANGE,  AND  ERROR 


A.  INTRODUCTION 

In  this  section,  we  consider  the  inference  of  the  concentration  and 
temperature  from  fan  beam  projection  data.  As  introduced  in  Section  I,  the 
projection  data  are  defined  by 


P(r,6)  =  -  In  [I(r,B)/I  ]  =  -  InR 

and  are  related  to  the  flow  cross-sectional  map  F(x,y) 

L 

P(r,6)  =  J  F(x,y)d£ 


where 

F(x,y)  =  q(vQ)nJ(x,y) 


(78) 


(79) 

(80) 


where  q(v)  is  the  absorption  cross  section,  nj  is  the  number  of  molecules 
in  rotational  level  J  of  the  molecule  under  investigation,  and  R  is  the 
ratio  of  intensities.  The  concentrations  nQ  and  nj  are  related  via  the 
Boltzmann  factor  when  the  gas  is  in  thermodynamic  equilibrium 


gjexp(-E/kT) 
nJ  =  no  Q(T} 


(81) 


where  gj  is  the  statistical  weight  of  energy  level  Ej  and  Q(T)  is  the  par¬ 
tition  function.  The  rotational  energy  level  diagram  of  a  typical  diatomic 
molecule  (N2,  CO)  is  shown  in  Fig.  24.  To  the  left  of  the  energy  levels  is 
a  scaled  plot  of  the  Boltzmann  factor  [Eq.  (81)],  indicating  the  relative 
distribution  of  nQ  among  the  various  levels  for  three  temperatures. 

If  the  flow  is  sampled  at  two  different  optical  frequencies 
v.j  and  Vg,  with  corresponding  absorptions  originating  from  levels  and 
J2,  with  cross  sections  q^  and  q2,  then  two  field  functions  F^(x,y)  and 
F2(x,y)  can  be  obtained  by  reconstruction  from  the  projection  data 
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Fig.  24  Energy  Level  Diagram  of  lZ  Ground  Electronic  and  Vibrational 
State  of  a  Typical  Diatomic  Molecule.  The  relative  population  in 
the  rotational  levels  are  shown  by  the  Boltzmann  curves  to  the 
left  for  three  temperatures. 


P^v-j)  and  P2(v2).  Substituting  Eq.  (81)  into  Eq.  (80),  writing  this 
expression  separately  for  each  absorption  at  and  v2>  we  can  solve  the 
resulting  two  equations  simultaneously  to  obtain  expressions  for  tempera¬ 
ture  and  concentration  maps 


F1(x,y)Q(T) 

nU.y)  =  -a  T'  )a -  exp  [E  /kT(x,y)  ]  (83) 

°  qrvrgJ1 

Equations  (82)  and  (83)  form  the  basis  of  our  analysis. 

B.  LINEARITY  AND  DYNAMIC  RANGE 

We  first  consider  the  dynamic  range  over  which  measurements  of  tem¬ 
perature  and  concentration  can  be  made,  given  a  linear  range  over  which  the 
ratio  R  can  be  measured.  The  linear  range  of  the  detector  system  deter¬ 
mines  the  linear  range  of  R.  The  detector  system  consists  of  a  fast  photo¬ 
diode,  the  sampling  electronics,  and  the  analaog-to-digital  (A/D)  con¬ 
verter.  Photodiodes  operating  as  photoconductors  are  linear  over  a  minimum 
of  seven  decades.  Fast  sampling  electronics  can  be  obtained  that  are  lin¬ 
ear  over  four  decades.  Thus,  a  10  or  12  bit  A/D  will  limit  the  linear 
dynamic  range  of  R  to  10^  or  4  x  10^,  respectively. 

To  find  the  dynamic  range  of  nQ(x,y)  for  a  given  dynamic  range  in  R, 
we  note 

R2  3n 

Ano  -  l  if  If  dR 
R1 
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From  Eqs.  (78)  and  (79) 


J  F  da  =  -  InR 


Taking  the  derivative  of  this  equation  with  respect  to  R  and  interchanging 
the  order  of  integration  and  differentiation  on  the  left  side,  we  have 


or 


r  3F  , 

J  3R  da  = 


3F  _d 
=  '  da  ^Rj 


3R 


Noting 


_  1_  dR 
R2  da 


dR  _  d_ 

da  "  da 


-njqa 

e 


"  njOR 


we  find 


3F 

3R  " 


-njq 

R 


(85) 


Taking  the  derivative  of  Eq.  (83)  with  respect  to  F,  substituting  this 
result  and  Eq.  (85)  into  Eq.  (84),  and  using  Eq.  (81)  for  nj,  we  find,  on 
integration 


An(x,y) 

n0(*»y) 


-  In  (Rg/R^ 


=  8.3  (12  bit  resolution) 


(86) 


68 


which  is  the  fractional  dynamic  range  for  a  given  pixel.  Equation  (86)  is 
the  same  result  one  obtains  for  simple  absorption  experiments.  This  dy¬ 
namic  range  can  be  extended  by  an  order  of  magnitude  by  viewing  two  dif¬ 
ferent  spectral  lines  simultaneously.  It  is  not  difficult  to  find  two 
absorption  lines  with  vastly  different  peak  absorptions  within  a  relatively 
narrow  spectral  range  [see  transitions  listed  for  ( 1 . 5 )  in  Table  4,  for 
example] . 

Similarly,  the  dynamic  range  in  temperature  can  be  estimated  from 


R„ 


r  Tb  JT_  Ab  3T  3F2 

“  1  3F,  '  lT  8F2  8R2  2 


Up 


Using  Eqs.  (80),  (82),  and  (85),  we  find 


(87) 


(88) 


C.  ERROR 

In  this  subsection,  we  consider  noise  in  the  reconstructed  image  and 
the  propagation  of  this  error  to  the  concentration  and  temperature.  Tom¬ 
ography  has  a  considerably  higher  sensitivity  than  single  ray  absorption 
due  to  the  averaging  effect  of  the  many  views  (projections).  Bennett  and 
Byer  (Ref.  43)  have  studied  the  noise  in  the  reconstructed  image  for  fan 
beam  geometry  using  a  convolution  back  projection  algorithm  that  utilized 
the  sampled  Ramachandran-Lakshminarayanan  (RL)  kernel.  When  the  noise  in 
each  projection  is  uncorrelated  with  the  noise  in  any  other  projection  (but 
allowing  for  correlation  within  a  given  projection),  the  reconstructed 
noise  variance  is 
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w(p) 


(89) 


where 


2  Q 

n  -  3 


l8Nfa2R2 


w(p) 


2/o 

p/2 


(1  - 


»¥ 


(90) 


where  op  is  the  variance  of  the  projection  data,  p  is  the  distance  from  the 
center  of  the  source-detector  ring  to  the  reconstructed  pixel  divided  by 
the  radius  of  the  source -detector  ring  Rp,  and  a  is  the  angular  separation 
between  the  detectors.  Hence,  the  uncertainty  in  a  given  pixel  in  the 
reconstructed  field  is  less  than  that  in  the  projection  data  by  the  square 
root  of  the  bracketed  term  in  Eq.  (89).  Note  that  the  noise  is  geometry 
dependent  through  the  term  w(p).  Thus,  keeping  F  (=  RD/r)  as  large  as 
possible  minimizes  the  "geometrical"  noise. 

Reference  43  defines  the  equivalent  noise  absorption  across  a  pixel 
as  e  min,  where 


mxn 


(rms  noise  amplitude)  x  (normalized  resolution  for  the  RL  kernel) 


op  x  a 


(91) 


Now,  consider  the  two-frequency  case,  where  both  temperature  and  con- 

2  2 

centration  are  measured  simultaneously.  Given  the  variances  a  and  o  of 

J  u  v 

the  two  measured  parameters  u  and  v,  the  variance  ax  of  the  desired 
parameter  x  can  be  estimated  using  (Ref.  44) 


2(  3X-1 

TO 


2r3xi 


(—1 

U  v3u'  '  °v  '•av' 


(92) 


where  ouv  is  the  covariance.  If  we  assume  that  the  fluctuations  in  R^  and 
Rg  are  uncorrelated,  the  covariance  will  average  to  zero  and  the  tempera¬ 
ture  variance  is 
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(93) 


Calculating  the  derivatives  from  Eq.  (82),  we  find 


(94) 


(95) 


assuming  that  opT  is  no  greater  than  op  oT.  Calculating  the  derivative  of 
nQ,  we  find  1 


exp(Ej/kT) 


2-2  -  ir-.-x 

°F  Q  +  °T  F1  laT 


2  -2/3Q  QEJr 


kT 


+  a^QFJ^  -  -J’ 


1/2 


T  F  11 3T 


kT 


(96) 


D.  EXAMPLE 

In  this  subsection,  the  dynamic  range  and  error  expressions  are  eval¬ 
uated  for  the  OH  molecule  and  the  geometry  discussed  in  the  previous  sub¬ 
sections.  The  hydroxyl  radical  is  an  important  minor  species  in  hydro¬ 
carbon  combustion.  Absolute  transition  probabilities  are  available  for  all 
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2  2 

electronic  transitions  in  the  (0,0)  band  of  the  X  i/2  “  ^  z  system 

(Ref.  45).  The  absorption  lines  from  this  band  appear  near  3100  A,  a 
spectral  region  accessible  by  a  doubled  dye  laser.  In  this  spectral 
region,  ultraviolet  (UV)  fiber  optics  can  transmit  enough  power  to  saturate 
the  detection  system.  Within  the  power  handling  capability  of  the  fibers, 
the  range  of  concentrations  that  can  be  measured  for  a  given  species  can  be 
changed  by  changing  the  laser  power  into  the  instrument  and  adjusting  an 
attenuating  filter  in  the  IQ  reference  channel  to  keep  the  range 
centered.  Thus,  the  12  bit  dynamic  range  in  R  can  be  maintained  for  values 
of  R  that  may  cover  several  orders  of  magnitude.  The  range  of  con¬ 
centration  can  also  be  adjusted  by  the  choice  of  spectral  line.  The 
absorption  cross  sections  of  most  molecules  vary  one  to  two  orders  of 
magnitude  within  the  spectral  range  of  a  single  laser  dye.  This  is  cer¬ 
tainly  true  for  OH. 

The  structure  of  the  ground  state  of  OH  is  a  spin  split  doublet  with  a 
splitting  constant  of  Ay  =  139.21  -  0.275v  (Ref.  46).  The  structure  is  not 
well  described  as  either  Hunds  case  (a)  or  case  (b)  coupling,  and  inter¬ 
mediate  coupling  formulas  should  be  used  to  calculate  the  structure.  The 

2  2 

splitting  is  inverted  with  the  n1/2  (F2)  state  lying  above  the  n^/2  (F^ ) 

state.  There  are  12  rotational  branches  of  electric-dipole-allowed  transi- 

? 

tions  between  the  ground  state  and  the  A  Z  state. 

The  energies  of  the  first  10  rotational  levels  of  both  spin  states  are 
given  in  Table  4.  The  fractional  population  of  each  of  these  levels  is 
calculated  at  a  temperature  of  1000  K  using  Eq.  (81).  The  partition 
function  for  this  molecule  can  be  approximated  by  (Refs.  46,  47) 


Q(T)  =  (2S  +  1)(2  -  fin  )(1  +  e“hCu)e/kT)(— J-) 


=  148  at  1000  K  for  OH 


where  S  =  1/2  and  A  =  1.  The  derivative  of  Q(T)  is  needed  for  Eq.  (96) 
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Table  4.  Selected  Properties  of  OH3 


2 

ni/2 

(f2) 

J" 

F2 ( cm- 1 ) 

f(10-2)* 

AJt  j.^sec"1 ) 

x  (A) 

qeff(cm2) 

0.5 

126.5 

1.13 

6.0  x  10~16 

1.5 

187.5 

2.06 

°P12(1.5) 

7.36  x  104 

3106 

2.5 

288.8 

2.68 

^(1.5) 

3.85  x  105 

3100 

3.1  X  10“15 

3.5 

429.4 

2.92 

4.5 

608 

2.82 

5.5 

824 

2.48 

°Q22(5.5) 

6.25  x  105 

3093 

2.5  x  10-15 

6.5 

1078 

2.01 

^22(5.5) 

4.6  x  104 

3093 

1.6  x  10“16 

■  7.5 

1368 

1.52 

Pq  (5.5) 

«12 

7.99  x  104 

3112 

3.3  x  10‘16 

8.5 

1693 

1.07 

9.5 

2054 

0.704 

10.5 

2453 

0.436 

^ — •» 

in 

d 

T— 

•w* 

C\J 

CM 

6.22  x  105 

3106 

2.5  x  10~15 

2 

n3/2 

(F,) 

J" 

F-j  (cm-1 ) 

f(10~2)b 

ANAJJtJ„(J”) 

AJ(  j.^sec-1 ) 

x(A) 

%fr(om2) 

■  1.5 

0.0 

2.71 

PP1l(1.5) 

8.6  x  i05 

3082 

6.9  x  10"15 

2.5 

83.9 

3.60 

RQ21 (1.5) 

1.8  x  105 

3072 

7.2  x  10'16 

3.5 

201.9 

4.05 

RRn(1.5) 

8.1  x  104 

3072 

2.1  x  10'16 

4.5 

355.1 

4.06 

5.5 

543.0 

3.72 

6.5 

767.5 

3.14 

•  7.5 

1029 

2.46 

QP21 (7.5) 

4.66  x  104 

3090 

2.1  x  10"16 

8.5 

1324 

1.81 

RRn(7.5) 

2.6  x  105 

3064 

9.2  x  10*16 

9.5 

1655 

1.25 

10.5 

2019 

0.81 

^(10.5) 

6.2  x  105 

3099 

2.5  x  1015 

a 


Energies  and  A  values  from  Reference  45. 


Fractional  populations  calculated  for  a  temperature  of  1000  K. 
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.-/niX  (2S  +  1 ) (2  -  6  )k  hew 

=  - h^B - —  l1  +  exp(-hcwe/kT)(l  +  -j^)]  (98) 


=  0.15  at  1000  K  for  OH 


The  absorption  cross  sections  in  Table  4  were  calculated  from  the 
Einstein  A  coefficients  using  (Ref.  49) 


8  j  2 

I  ■)<v)dv  =  (r1)  k  #j„j. 


To  evaluate  Eq.  (99)  for  a  value  of  the  cross  section,  we  must  assume  a 
form  for  the  lineshape.  If  qQ  is  the  cross  section  at  the  line  center  v  , 
and  G(v)  is  the  lineshape  function,  then 


q(v)  =  qQG(v) 

For  a  Doppler  (inhomogeneously)  broadened  line,  G(v)  is  Gaussian,  i.e. 

Gq(v)  =  exp[(v  -  vo)2/y^] 

where  GD(vQ)  =  1,  jGD(v)dv  =  ^  ,  and  yG  is  the  e-1  halfwidth.  For  a  pres¬ 
sure  (homogeneously)  broadened  line,  G(v)  is  Lorentzian,  i.e. 

Gd(v)  =  Yp[(v  -  vQ)2  +  Yp]"1 


where  Gp(vQ)  =  1,  J  Gp(v)dv  =  ttyp>  and  yp  is  the  half  width  at  half  maxi¬ 
mum  (HWHM).  Note  that  yd  =  Yp//£n2.  An  OH  molecule  in  a  combustion  envir¬ 
onment  has  absorption  lines  that  are  primarily  pressure  broadened,  so 


For  laser  radiation  with  a  finite  bandwidth  centered  at  vq,  the  peak  cross 
section  is  not  realized  for  all  the  laser  radiation,  and  an  effective  (so- 
called  midband)  cross  section  equal  to  0.5  q0  is  used.  The  values  of  this 


cross  section  are  listed  in  Table  4.  The  values  for  A  in  Table  4  were 
taken  from  Reference  45.  The  linewidths  are  functions  of  T,  J,  and  the 
collision  environment.  An  average  value  of  0.05  cm“^  was  adopted  for  the 
halfwidth  (Ref.  50). 

We  consider  first  an  experiment  using  a  single  absorption  line.  For 
this  situation,  only  nj(T)  can  be  measured  for  each  pixel,  and  the  relative 
number  density  cannot  be  measured  unless  either  T(x,y)  =  constant  or  nj  is 
not  a  function  of  temperature.  The  first  case  is  physically  uninterest¬ 
ing.  Fortunately,  the  second  case  is  possible  because  most  molecules  (that 
are  amenable  to  this  technique)  have  a  level  from  which  a  spectrally  acces¬ 
sible  line  originates  whose  population  is  relatively  insensitive  to  changes 
in  T.  For  example,  the  absorption  from  any  line  originating  from  the  Ej  = 

0  energy  level  is  a  function  of  temperature  only  through  the  partition 
function.  The  partition  function  to  first  approximation  is  a  linear  func¬ 
tion  of  T  [Eq.  (97)],  and  a  200  K  variation  in  T  introduces  a  20%  uncer¬ 
tainty  in  the  inferred  values  of  the  concentration  for  temperatures  near 
1000  K.  In  Table  4,  the  PP^(1.5)  and  RR^(1.5)  both  originate  from  this 
level  and  have  absorption  cross  sections  that  differ  by  nearly  two  orders 
of  magnitude,  allowing  both  high  and  low  concentration  environments  to  be 
studied.  The  dynamic  range,  however,  is  still  limited  by  Eq.  (86).  For  a 
given  temperature  range,  higher  accuracy  may  be  attainable  by  choosing  a 
rotational  level  where  the  variation  in  the  partition  function  due  to  tem¬ 
perature  is  offset  by  the  change  resulting  from  the  exponential  function 
[Eq.  (83)]. 

Using  Eq.  (91),  the  noise  level  in  the  density  measurement  can  be 
estimated.  If  we  assume  that  the  ratio  of  intensities  R  (=  I / IQ )  can  be 
measured  with  10  bit  accuracy,  the  relative  error 


is  0.1%.  From  Eqs.  (78)  and  (92),  the  projection  variance  is 
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(102) 


3P 

°P  =  °R  3R  "  R  R 

The  projection  error  is  related  to  the  error  in  the  field  function  error 
through  Eq.  (89).  We  consider  our  optimum  geometry  for  a  4  in.  diameter 
flow  (r  =  5  cm)  to  have  n^  =  96,  F  =  3.86,  and  choose  an  average  r/2  = 
2.5  cm  to  obtain  a  value  for  p  of  0.1275.  Then,  from  Eq.  (90),  w  =  1.01 
and  Op  =  0.31  op.  Allowing 

6  .  =  n  .  qa,  =  n  qFLa 

mm  noise^  res  noise  D 


we  have,  from  Eq.  (93) 


n 


noise 


qRr 


(103) 


For  the  (1.5)  line,  nnoige  =  2  x  109  cm-3  =  7  x  10  11  Amagat.  For  a 

4  in.  flow  diameter,  the  maximum  detectable  density  is  10 cm  J  for  an 
attenuation  of  ~  101*.  Because  of  limited  dynamic  range,  the  minimum 
detectable  intensity  is  -  1013.  Thus,  this  noise  level  is  ~  0.02%  of  this 
lower  level,  a  value  considerably  below  other  sources  of  noise,  such  as 
laser  speckle. 


The  more  important  experiment  is  the  probing  of  two  different  spectral 
lines  of  a  single  species.  By  the  proper  choice  of  lines,  the  dynamic 
range  of  concentration  can  be  increased  by  nearly  an  order  of  magnitude, 
and  a  spatial  map  of  the  temperature  also  can  be  recovered.  To  increase 
the  dynamic  range  of  nQ,  two  lines  with  an  absorption  cross  section  that 
differ  by  a  factor  of  ~  10  should  be  chosen.  According  to  Eq.  (88),  the 
dynamic  range  of  T  varies  directly  with  the  energy  difference  between  the 
lower  levels  of  the  two  lines.  From  Eq.  (9*0,  this  is  also  true  for  o^.. 
Because  the  dynamic  range  is  large  for  T,  closely  spaced  levels  should  be 
selected.  The  dynamic  range  of  T  can  be  increased  by  choosing  lines  that 
originate  from  opposite  sides  of  the  peak  of  the  Boltzmann  curve  so  that 
the  two  logarithmic  terms  in  the  last  bracket  of  Eq.  (88)  add.  Finally, 
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note  that  the  number  density  can  be  calculated  from  either  or  R2.  The 
variance  of  nQ  is  decreased  by  choosing  an  absorption  from  an  energy  with  a 
low  value  of  exp(Ej/kT)/gj  and  with  a  large  absorption  cross  section.  With 
these  criteria  in  mind,  we  select  the  PP22  (1.5)  and  PQ.j2  (5.5)  lines. 

Equation  (88)  is  evaluated  to  obtain  the  dynamic  range  in  T  assuming 
that  the  range  of  the  ratios  for  the  J2  line  may  take  the  full  12  bit 
range,  but  the  range  of  ratios  for  is  a  factor  of  10  less.  Accordingly, 
we  find  AT  =  7850  K.  The  dynamic  range  of  nQ  is  ln(4  x  10^)  x  {q[PP22(1 .5)]/ 

D 

q [  Q 12(5-5)3}  »  78.  The  variance  in  temperature  is  given  by  Eq.  (9*0. 

Noting  F  =  qn,  and  evaluating  for  a  partial  density  of  10"^  Amagat,  we  find 
Orp  =  7.6  for  eR  =  0.1/6.  Finally,  evaluating  Eq.  (92)  for  the  above  condi¬ 
tions,  we  find  oR  =  2  x  lO1^  =  8  x  10“^  Amagat — an  accuracy  of  0.8/L 
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